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Abstract. The relation between the dynamical regimes (weak and strong coupling) 
and entanglement for a dissipative quantum - dot microcavity system is studied. 
In the framework of a phenomenological temperature model an analysis in both, 
temporal (population dynamics) and frequency domain (photoluminescence) is carried 
out in order to identify the associated dynamical behavior. The Wigner function and 
concurrence are employed to quantify the entanglement in each regime. We find that 
sudden death of entanglement is a typical characteristic of the strong coupling regime. 
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1. Introduction 

In the last few years, the study of microsystems that involve the interaction between 
an active medium and confined light, has made possible the observation of interesting 
phenomena in two dynamical regimes: weak and strong coupling [1, 2, 3]. In the first 
regime, spontaneous emission control was successfully realized (Purcell effect) [4, 5]; in 
the second one, several groups are now searching for coherent polaritonic phenomena 
such as lasing [6] or condensation [7, 9], which can open applications in quantum infor- 
mation and quantum optics [8, 10, 11, 12]. In addition, some recent theoretical works 
have shown that the dynamical properties for a coupled quantum dot - cavity system, 
may be described using a simple dissipative model [13, 14]. 

Our purpose in this work is to study the relations between weak and strong coupling 
regimes with the dynamical exciton - photon entanglement, by using the concurrence 
measure and the Wigner quasiprobability function. Despite the fact that the Wigner 
function depends only on the photonic state, it has been shown that it can be used as 

a qualitative criterion for determining the separability of the quantum exciton - photon 
state [15]. In reference [16] a similar study has been carried out but studying 2 quantum 
dots and only in the stationary limit. 

The paper has been written as follows: section 2 contains the theoretical frame- 
work supporting our model, in section 2.1 we describe the system and the dissipative 
processes that models its dynamics using a master equation. In section 2.2 we explain 
how the photoluminescence spectrum is obtained by using the quantum regression the- 
orem. In section 2.3, using a simple phenomenological model we include temperature 
effects in the quantum dot gap. In section 2.4 we review the concurrence and Wigner 
function concepts, and their connection with the photon - exciton entanglement. Hence, 
the dynamical regimes arc characterized employing the numerical integration results of 
the master equation for two different cutoff conditions in the photon number, in both 
cases, showing a good agreement with experimental data. Once the regime is charac- 
terized, a dynamical description of the entanglement is analyzed and we show collapses 
and revivals of this quantity as a function of the dissipative parameters involved in the 
model. Furthermore, we establish the usefulness of the Wigner function criterion to 
detect separability in a multi-state system where the concurrence criterion can not be 
used. Finally, some conclusions are given in the last section. 

2. Theoretical background 

We are interested in the evolution of a quantum dot interacting with a confined mode of 
the electromagnetic field inside a semiconductor microcavity. In these systems, quantum 
states associated to the matter excitations, the so-called excitons, are bound states 
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resultant of the Coulomb interaction between electrons in the conduction band and 
holes in the valence band. This quasi-particle exhibits a complete discrete excitation 
spectrum, however, in this work we will only consider the first two levels of this set, the 
ground \G) (no excitation, i.e, electron in the valence band) and excited |X) (exciton) 
states. This assumption is based upon the fact that the ground and first excited states 
of a multilevel model involving Coulomb interaction, are mainly filled in the dynamical 
evolution of the system [17]. A possible experimental realization of this model, can be 
implemented in a pumped system with polarized light [18] and slow spin flip mechanisms. 
The photonic component will be treated as a single electromagnetic mode. The vahdity 
of this usual assumption is subjected to the existence of well separated modes in energy 
inside the cavity [19]. The last condition amounts to say that for example the radii 
of the micropillar is small since the energy separation of the modes increases when the 
radii is decreased. We will employ a Fock state basis |n), which ought to be truncated 
in order to implement computationally the dynamics. 

Light-matter interaction is described by the Jaynes-Cummings Hamiltonian, 



where a — and cr^ = ^^e exciton ladder operators and a (a^) is 

the annihilation (creation) operator for photons, cu and a;o are the exciton and photon 

energy, respectively, and g is the coupling constant and we have set h—1. 
Wc also define the detuning between the exciton and photon frequency as A = u — ujq. 
Under these considerations the system is Hamiltonian and completely intcgrablc [21]. 
However, real physical systems are far away from this simple description; dissipative 
effects play an important role in the evolution of the system. Indeed, if no losses were 
considered in the system, no measurements could be done since light would remain 
always inside the microcavity. 

2.1. Dynamics 

The whole system-reservoir hamiltonian can be splitted in three parts. One is for the 
system we are considering, namely, the photons of the cavity and the exciton. The second 
one is the hamiltonian of the reservoirs, which is made of electron-hole pairs, photons 
and phonons and finally a third part which is a bilinear coupling between the system 
and the reservoirs. The explicit form of the system-reservoir interaction for this model 
can be found in [13]. After tracing over the external reservoirs degrees of freedom and 
assuming the validity of the Born-Markov approximation, which requires weak coupling 
between the system (cxciton+cavity photons) and the reservoirs one arrives to a master 
equation. It has recently been found [20] that non-markovian dynamics is relevant for 
high pumping intensities. 

The master equation we have found accounts for three different processes, namely, 
coherent emission (k), external pumping (P) and spontaneous emission (7). The master 
equation we shall consider is [13, 22, 23]: 




(1) 
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By changing the values of the free parameters in this model, two different dynamical 
regimes are reached: weak and strong couphng. Possible transitions between these two 
regimes, can be achieved as loss and pump rates are modified. The dynamical behavior 
of the system as well as the size of the basis employed are governed by the competition 
of the time scales involved in the model. The time scales associated with the non 
conservative processes included are given by: rp = l/P, = 1/k and = I/7, 
whereas, the interaction time scale is Tg = 1/g. We address now to the interpretation 
of these time scales. In figure 1 three possible cases arc shown whose dynamics can be 
described using a basis including up to one photon (cases (A) and (B)) and more than 
one photon (case (C)). 

In case (A) , for instance, the relation k> P ^ g holds, and the system is operating 
in weak coupling regime. When an exciton is pumped during the typical time scale Tp, 
the elapsed time until it recombines into a photon is given by the (largest) scale Tg. Bear 
in mind that throughout this period no further excitation can be done over the quantum 
dot, because the Pauli exclusion principle does not allow an additional excitation. The 
photon generated in this way, quickly leaves the cavity due to the small time scale 
in which it can inhabit the cavity. The latter mechanism applies for any photon living 
the cavity. 

Hence we see that there is no chance for Rabi oscillations in the dynamical evo- 
lution. If the initial condition for the electromagnetic field is the vacuum state, the 
photon mean number is never expected to be greater than one. 

A similar reasoning can be done in the case (B), where g > k ^ P. In this case each 
photon created by the interaction can be re-absorbed and produce Rabi oscillations 
(i.e., strong coupling), before it leaves the cavity. A basis inchiding up to one photon 
is enough to describe the system because the characteristic time related to losses is 
by much smaller than the pump one. The important issue here is that the dynamical 
situation is clearly different to the previous one. 

Finally, in case (C), where g > P ^ k, we have a situation in which the photons can be 
efficiently stored within the cavity. Indeed, since the coherent emission rate is small, its 
characteristic time scale is large enough compared with the associated excitonic pump 
and interaction rates. Photons created by exciton recombination remain in the cavity 
a long time before they leak through the cavity mirrors. This case is a clear example 
where the multiphoton basis must be implemented. 
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(A) ^ 1 > K>P:>g 
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(B) ^ \ > g>i^:^P 

Tg Tk Tp 

(C) ^ 1 > g>P:^K 

TgTp Tk 

Figure 1. Schematic diagraik of some relations among the typical time scales 
associated with the effective rates involved in the master equation. Cases (A) and 
(B) show situations in which the basis containing up to one photon is enough for 
describing the dynamics. On the other hand, case (C) shows a set of parameters where 
this description fails, and the requirement for a larger basis arises. See the text for 
additional details. 

2.2. Photoluminescence 

The Fourier transform of the first order correlation function directly gives the 
photoluminescence spectrum of the system [21], 

/oo 
{a^t + T)a{t))e^^dT. (3) 
-oo 

Note that in this expression a knowledge of the time-dependent expected value for 
the product of creation and annihilation operators, is needed in order to compute the 
photoluminescence spectrum. However, non analytical expression for such expectation 
value is available for our system. This problem can be solved by representing the 
operators a (a^) in the interaction picture, and then using the quantum regression 
theorem [22]. This theorem states that given a set of operators Oj satisfying, 

^mt+T)) = j2Ljkmt+T)), (4) 

k 

then 

^{Oj{t + T)o{t)) = Yl Li^mt + T)o{t)). (5) 

k 

for any operator O. The validity of this theorem holds whenever a closed set of 
operators are involved in the dynamics in the Markovian approximation. Unfortunately, 
representing creation and annihilation operators in the interaction picture does not lead 
to a complete set. It is necessary to add two new operators in order to close the system. 
The final set of equations for the operators are: 
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a^^,^Xt) = |Xn + l)(Xn|e'('"-^^* 
ai{t) = \Xn){Gn\e''^^ 

= \Gn + l){Xn - l|e*(-2^)*. (6) 

In the framework of the previous discussion regarding the dynamical time scales, we 
can consider the simplest model involving all dissipative and interaction processes, that 
is, we shall study a problem having just three levels: ICO), \X0) and \G1). Taking 
into account this consideration we only need to define two operators whose dynamical 
equations are: 



^(4o(^)) = - (f + ^ - ^(^ - A)) {aUt}) + wi^^M 

j^oit)) = ^9{aUi)) - - (^iW), (7) 

thus by using the quantum regression theorem we can write the dynamics for the delayed 
operators as, 

X^±( (4o(^ + ^)) 
dr\ {alitor)) 



;| + P) + ,(^_A) ig \ I {al,{t + T)) 
ig -f±l+ia; M {al{t + T)) 



AX. 



The last linear system has the following formal solution 

X{t + r) - e^^X{t) = Be^^B-'X{t), (8) 

where A — BAB"^ and A = diag{A+, A_} is a diagonal matrix containing the eigenvalues 
of A. 

The eigenvalues A± are directly related with the peaks u!± and widths r± of the 
spectrum, 

u;± + ir± = iX±. (9) 



2.3. Temperature effects 

In order to make a more realistic description of the experimental data we also include 
a temperature dependence in the model. It is important to explicitly point out that 
the developed master equation describes the system at zero temperature, that is, the 
states of the reservoirs considered in the derivation of the master equation (2) are zero 
temperature states. This assumption can be justified as follows. 

First notice that when a finite temperature reservoir of photons is considered 
in the derivation of the master equation one obtains two Lindblad terms [24]. One 
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that accounts for thermally induced absorption and that is proportional to the average 
number of photons of the reservoir and another term that is proportional to the average 
number of photons of the reservoir plus one and that accounts for spontaneous emission. 

The average number of photons of a thermal reservoir at temperature T is 
l^^^ojo/kBT _ j^-j rjj^g standard experimental values for the photon energy in a mi- 
cropillar are ^ 1 meV whereas for temperatures of the order of 10^ K the thermal 
energy ksT ~ 1 meV and thus the average number of photons is N{u) ~ e~^°^ ~ lO"^"*^. 
We see that effects due to finite temperature in the master equation for the system con- 
sidered here are quite small. 

Summarizing, the above discussion, we can add thermal effects on the system 
through its effective parameters, such as the quantum dot energy gap and the cavity 
refractive index, without considering the fundamental issues introduced into the model 
by the dephasing processes. First, let us consider the temperature photon energy 
dependence. Usual microcavities are built from GaAs/AlGaAs layers, for these cavities 
the resonant wavelength depends on the refractive index as A = Xaii/n, where Aair is the 
light wavelength in vacuum. We can then modify the resonant frequency by changing 
the refractive index n{T) [12]. Experiments have shown that this index has an almost 
linear temperature dependence within the range from K up to few hundred Kelvin. In 
reference [25] has been shown that the refractive index can be modeled with the simple 
formulae: 



Where a ~ 10~^ K~^. It is seen that the corrections to the wavelength due to the 
temperature are rather small, and because of this will not be considered in this work. 
On the other hand, the temperature effects, related to the active medium can be included 
through the energy gap in the quantum dot. We shall use here the Varshni model [26], 
which fits the band gap thermal dependence in the low temperature region. A more 
detailed discussion about the validity of this model, can be found in reference [27], 



In this work we will use two different sets of experimental data for InGaAs quantum 
dots, for which the Varshni parameters were fitted. 

2.4- Concurrence and Wigner Function 

Determining entanglement between two quantum mechanical systems is a complicated 
task when the systems involved have several degrees of freedom, that is, when the 
basis representing the density matrix has more than two pairs of states. However, if 
each system is completely described with two levels (i.e. each system is a qubit) the 



n{T)^no{l + aT). 



(10) 



uj{T) = Eg{0) 



(11) 



Characterization of dynamical regimes and entanglement sudden death. 



8 



situation becomes easier; indeed, the entanglement can be obtained from the concurrence 
expression [28], 

C(p) = max{0, Ai - A2 - A3 - A4}, 

where {Aj} are the square root of the decreasingly ordered eigenvalues of the 
positive-definite non-Hermitian matrix pp with, 

p= {ay^ay)p*{ay^ay), 

and p* is the complex conjugated matrix representation of p in some suitable basis. 

For instance, note that if our system is described by using the basis {\G), \X)} 
{|0), |1)}, (where |0), |1) are photon Fock states), the exciton - photonic field system 
can be thought as two interacting qubits. Once a (numerical) solution of the master 
equation is obtained, it is straightforward to compute the concurrence following the 
previous recipe. It is to be noticed that the validity of this basis is constrained to 
the dynamical regime where this cutoff holds. When the basis is larger than the one 
considered earlier, the concurrence measure is not applicable and hence a new criterion 
must be established. Despite the Wigner function only depends on the photon state, 
it has been demonstrated that it yields qualitative information about the separability 
between the exciton and photon parts of the global quantum state [15]. The Wigner 
function for the photonic field can be easily computed as [29], 

W(a) = 2 Trp[D(-a)p^^^D(a)Pf], (12) 

where Pj = c''^"^" is the field parity operator and D{a) = e""^""**^ is the displacement 
operator, p^^^ is the density operator of the photons and Trp is the trace operation 
in Fock space. Our system involves both excitonic and photonic states, hence this 
definition is not directly applicable, and a previous differentiation among the possible 
combinations has to be done, that is, we have to separately consider the Wigner matrix 
elements Wij [15], 

Wijia) = 2 TTp[D{-a){i\p\j)Dia)Pj], 

where indexes i.j run over excitonic states IG)}, and p is the density operator of 

the whole system. Notice that {i\p\j) is an operator that acts only in the state space of 
the photons so that the operations in the above equation are well defined. 

Note that if the system is separable at a given time then the matrix elements Wxx 
and Wqg have the same shape in phase space, leading to a separable system, which is 
itself related to a vanishing concurrence situation. Indeed, recognizing that each Fock 
state has a well defined signature within the phase space, it is possible to identify the 
predominant photonic state by observing the Wigner function. A more detailed analysis 
shows that if the quantum state is separable, 

I 
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(13) 



where 1^(q;) = Yli^ai'^)- Therefore, if the system is separable at a given time all 
the matrix elements Wij of the Wigner matrix have the same shape for every i,j. In 
opposition to the quantitative concurrence measure, the Wigner function only shows 
qualitative information about the entanglement. 

3. Results 

3.1. Weak coupling regime 

In the weak couphng (WC) regime the relation [23] 



holds, hence by selecting a strong enough emission rate WC is obtained. We choose 
g = 15 /xeV, 7=1 ^eV, n = 85 fieV and P = 20 fieV corresponding to the case (A) 
of figure 1. As we stated before, this situation can be described with the reduced basis 
including up to one photon. 

We stress that relation (14 ) is derived under the assumption that there is no incoherent 
pumping over the system. The effects of the pumping have been studied in detail in 
reference [30]. 

We have taken a photon decay rate of k = 85 /leV and a cavity mode of fre- 
quency uj = 1296.11 meV, so that the quality factor of the cavity we are considering 
is Q ~ 15000. The spectra for this set of parameters is shown in figure 2. These val- 
ues were chosen following the results of reference [11] where strong coupling in a single 
quantum dot - microcavity system is reported. 

In figure 2 it is seen that at approximately T = 12i^ the peaks of the cavity and exciton 
coincide. The parameters used for the temperature dependence where -E'g(O) = 1299.6 
meV, a = 0.81 meV/K, /3 = 457.6 K. 

We pointed out that weak coupling regime could be reached by setting an emission 
rate higher than the coupling constant. Due to the small time the photons spend in the 
cavity the chance of them to interact with the exciton is quite small and because of this 
no oscillations appear in the populations as seen in Figure 2. 

3.2. Strong coupling 

The strong coupling (SC) regime is even more interesting since it enables the existence 
of polaritons. From an experimental point of view it is more complicated to obtain due 
to the fine scales of the variables involved. This regime is also studied with the same 



16g^ <{k- 7)2, 



(14) 
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Figure 2. (left panel) Crossing emission peaks in the WC regime computed with a 
simplified model (solid line), the whole model using a 40 photonic level basis (squares), 
(right panel) Time evolution of the populations and the average photon number in the 
WC regime, computed using the simphfied model, 3=15 /ieV, 7 = 1 /^eV, k = 85 ^eV 
and P = 20 /leV. The average photon number was computed using 40 Fock states. 



three level model described above and with a basis of 40 Fock states but now taking, 
= 35 /ieV, K = 25 /xeV, 7 = 1 /xeV, P = 1 /leV, Eg{0)= 1299.6 meV, a= 0.81 meV/K, 
/?= 457.6 K and a cavity resonant frequency of u = 1299.35 meV. 
It is clear that in this regime oscillations in the populations are observed since the time 
a photon can live in the cavity is long enough for it to interact with the exciton, i.e., 
1/g = Tg> T^ = 1/k. 




12 14 16 50 100 150 200 250 300 

T(K) " t (ps) 



Figure 3. (left panel) Emission peaks in the SC regime. Anticrossing can be observed 
at T = 12 K. Solid line was computed by using the simple three level model, meanwhile 
dots are the peaks obtained with a basis involving 40 Fock states, (right panel) Time 
evolution of populations and the average photon number for the strong coupling regime, 
again only the average photon number was computed with 40 Fock states. Differences 
with the weak coupling regime are clear. 

Peaks positions are shown in fig. 3. 

Note that for detunings far from the resonance, a small shift of the cavity peak 
appears, this effect is not observed in the WC regime. This shift is due to the strong 
interaction between matter and light. 
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The closest approximation of the two emission peaks (cavity and excitonic) is reached 
when T — 12 K, where the separation (Rabi sphtting) is 

R = 68.89 neV, 

Another important fact which enables us to conclude that we are dealing with a system 
in SC regime is that the following condition holds 

That is the Rabi splitting is approximately two times the coupling constant [23]. 
Finally to determine if the reduced basis was enough to describe the system we computed 
the average photon number (^40) using a basis of 40 Fock states. This is information 
is in the right panels of figures 2 and 3 and can be compared with the average photon 
number using the reduced basis: (rii) = pciGi- It is seen that both descriptions are in 
complete agreement. 

3.3. Concurrence 

In the previous section a study of the relation between dynamics and photoluminescence 
was made, leading to a direct identification of strong and weak coupling regimes. Now 
we want to study the behavior of the entanglement in these regimes. In order to do so 
we now turn our attention to the simplified system described in section 2.4. First we 
study the system in strong coupling regime with a set of parameters as those described 
in 3.1; however, since we are mainly interested in the effects that dissipation has on the 
entanglement we set A = 0. Since spontaneous emission processes are small as com- 
pared with the other effects, we also set 7 = 0. Taking g — 25 /leV we obtain the results 
in the upper panels of fig. 4. In the middle left panel the evolution of the concurrence 
as a function of detuning is shown. Notice that when the system moves from strong 
coupling to weak coupling the entanglement becomes larger, leading to a non vanishing 
value of the concurrence for all times. 

The range of parameters for this situation was carefully chosen since great varia- 
tions of them lead to unphysical results. Indeed, if k is increased the cavity eventually 
will get empty, if P is increased the basis is not large enough to describe the system, 
since the mean number of photons becomes larger than 1. 

Results on the evolution of concurrence show that it does not distinguish between 
the two dissipative processes (parameters) considered {k, and P). The reason for this 
behavior is that the concurrence of our system depends essentially of the difference 
between the absolute value of the coherence term pcixo and the square root of the 
product of the populations pcoGo and pxixi- Notice for example that in figure 4 
panel (e) the concurrence vanishes for precisely the times when |pgixo| < y/PaoGoPxixi 
and that the difference between the concurrence and |pgixo| in their minima grows 
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Figure 4. Concurrence time evolution as a function of the dissipative parameters: 
(a) changing k and fixing P = 3 /icV, A = , (b) changing P and fixing k = 3 /^eV, 
A = and (c) changing A and keeping P ~ 2 fieV, k = 3 fxeV . In (d) time evolution 
of linear entropy {M{p) = 1 — Tr(p^)) for the last set of parameters. In figures (e) and 
(f) we plot the concurrence (C(p)), linear entropy (M(p)), y^pcoGopxixi and |pxoGi| 
for A = eV and A = 0.1 eV respectively with k = 3 /icV, P — 2 /leV. Dynamics was 
solved with the initial condition pxoxo = 1 and the coupling constant g = 2b fieV in 
all cases. 
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Figure 5. Revivals time 5ti (left panel), 5t2 (right panel) in ps as functions of k and 
P. As dissipative factors increase the elapsed time in the concurrence revivals becomes 
longer. 



as the coherence y/pcoGoPxixi grows as seen in panel (f). In the dynamics of the 
density matrix we are considering it is seen that the coherences decay exponentially 
with rates proportional to k + P and that on the other hand P tends to increase the 
population of pxixi and n tends to increase the population of Pgogo so that the difference 
IpgixoI — y/PGOGoPxixi docs uot distinguish the two processes. 

Even more interesting is the fact that as the dissipation increases the zones where 
concurrence vanishes (the so called Entanglement Sudden Death [31]) become wider, 
revivals in the concurrence become more separated as dissipation increases. Note that 
when dissipation increases the maxima of the concurrence are less defined and eventually 
disappear. In order to quantify these two effects we compute the temporal length of 
the first two collapses of the concurrence {6ti is the time interval of the first collapse 
and 6t2 the time interval of the second one). The results obtained are plotted in 
fig. 5. It is clearly seen that as the non hamiltonian effects become important the 
length of both intervals gets longer and that for any k and P the second interval 
St2 is wider than the first one 6ti. This behavior can be understood as follows: as 
explained in the last paragraph the degree of entanglement between the subsystems 
is the difference |pgixo| — y/PcoGoPxixi and the dynamical behavior of the coherence 
term is an oscillatory function times an exponential decaying function with decay rate 
proportional to the sum of k and P. In order to have entanglement the absolute value of 
IpgixoI must be greater than ^/pcoGoPxixi- Now notice that as P or k are increased the 
oscillations gets more damped. Because of the damping the amplitude of the oscillations 
is smaller and the absolute value of the coherence has to be in a time t closer to a 
maxima in order to be greater than y/pcocoPxixi, that is, the finite time where there 
is no entanglement approaches to the time interval between to successive maxima of 
IpgixoI as k or P is increased. To understand the fact that the first interval of sudden 
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Figure 6. Time behavior for Wigner function matrix elements at different concurrence 
regimes. Notice that in the times of non vanishing concurrence {ti and ^3) the shape 
of the matrix elements is quite different, in particular the non diagonal elements can 
not be obtained by multiplying the trace (which is a real valued function) since they 
have more than one lobule. On the other hand in the vanishing concurrence case (^2) 
the shape is the same. 

To further quantify the separability of the quantum state of the whole system we 
can define the following quantity: dWij{a) = Wij{a) — 'vv(o')^ W{a) if the state is 
separable then 6Wij{a) = Va. Actually when one calculates 6Wij{a) in ti and 
SWij{a) is bounded as follows < |i5Wij(a)| < 3 and in t2 is has the following bound 
< < 0.04. The times ti, <2 and ^3 are indicated in figure 4. 

death Sti is smaller than the second one St2, simply notice that for the times in the 
second interval the amplitude of |pgixo| will be smaller than in the first one, so that the 
times when the entanglement will be different from zero will be closer to the maxima of 
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IpgixoI resulting in a wider interval. 

This behavior suggest that by manipulating experimentally accessible parameters 
such as the pumping rate and the quality factor of the cavity one can in principle make 
coherence control in the system. Finally, it is worth to mention that due to the coupling 
with external reservoirs the system evolves to a non pure state as evidenced in figure 
(4) where we plot the linear entropy. 

Now we compute the elements of the Wigner function Wgg and Wxx (see fig. 6) 
in representative zones of the concurrence function for typical values of the parameters: 
K — S fjLeV and P — 2 //eV, at the times ti, t2 and ^3 indicated in figure 4. Note that in 
the regions where concurrence goes to zero, phase space is similar (up to a multiplicative 
constant) as suggested by Wxx, Wgg, ^^{Wgx) and \m{WGx), so that the photonic 
states are similar for both ground and excited excitonic states and hence the photonic 
states are separable from the excitonic part leading to a vanishing entanglement. On 
the other hand for non vanishing concurrence it is clear that photonic states can not be 
separated. 

4. Conclusions 

We have built a phenomenological model that is able to describe both weak and 
strong coupling an that accounts for temperature effects in a microcavity quantum dot 
system. We have used concurrence criterion and the Wigner function to carry out 
the entanglement analysis related with the dynamical regimes in a simple dissipative 
model. The strong coupling regime shows a periodical disentanglement that depends 
on the dissipation rates. On the other hand, the weak coupling regime shows no 
complete dynamical losing of entanglement. This relation between dynamical regimes 
and entanglement sudden death has been investigated and we have shown that on the 
one hand both the coherent emission (k) and the incoherent pumping (P) affect the 
time windows where there is no entanglement. In the same way, as time goes by this 
windows become wider due to the damping in the dynamics of the non diagonal terms 
of the density matrix caused by n and P. 

Finally, we would like to point out that our results suggest that one can control the 
entanglement between the subsystems by manipulating external parameters such as the 
pumping rate P and the cavity quality factor (which is related to the emission rate k). 
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